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We use distorted wave electron scattering calculations to extract the weak charge form factor 
Fw{ctj, the weak charge radius Rw, and the point neutron radius _R„, of '^"^Pb from the PREX parity 
violating asymmetry measurement. The form factor is the Fourier transform of the weak charge 
density at the average momentum transfer q = 0.475 fm"\ We find Fw{q) = 0.204 ± 0.028(exp) ± 
0.001 (model). We use the Helm model to infer the weak radius from Fw{q)- We find Rw = 
5.826 ± 0.181(exp) ± 0.027(model) fm. Here the exp error includes PREX statistical and systematic 
errors, while the model error describes the uncertainty in Rw from uncertainties in the surface 
thickness a of the weak charge density. The weak radius is larger than the charge radius, implying 
a "weak charge skin" where the surface region is relatively enriched in weak charges compared 
to (electromagnetic) charges. We extract the point neutron radius _R„ = 5.751 ± 0.175 (exp) ± 
0.026(model) ± 0.005(strange) fm, from Rw- Here there is only a very small error (strange) from 
possible strange quark contributions. We find Rn to be slightly smaller than Rw because of the 
nucleon's size. Finally, we find a neutron skin thickness of Rn — Rp — 0.302 ± 0.175 (exp) ± 0.026 
(model) ± 0.005 (strange) fm, where Rp is the point proton radius. 

PACS numbers: 21.10.Gv, 25.30.Bf, 24.80.+y, 27.80. +w 



Parity violating elastic electron scattering provides a 
model independent probe of neutron densities, because 
the weak charge of a neutron is much larger than the 
weak charge of a proton [T . In Born approximation, the 
parity violating asymmetry Ap^ , the fractional difference 
in cross sections for positive and negative helicity elec- 
trons, is proportional to the weak form factor F^. This 
is very close to the Fourier transform of the neutron den- 
sity. Therefore the neutron density can be extracted from 
an electro- weak measurement [T]. However, one must in- 
clude the effects of Coulomb distortions, which have been 
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accurately calculated [1], if the charge density Pch [3] is 
well known. Many details of a practical parity violat- 
ing experiment to measure neutron densities, along with 
a number of theoretical corrections, were discussed in a 
long paper jj. 

Recently, the Lead Radius Experiment (PREX) mea- 
sured Apy for 1.06 GeV electrons, scattered by about 
five degrees from ^'^^Pb, and the neutron radius i?„ was 
extracted [5]. To do this, the experimental Ap^ was com- 
pared to a least squares fit of i?„ as a function of Apy , pre- 
dicted by seven mean field models [6j (see also In the 
present paper, we provide a second, more detailed, anal- 
ysis of the measured Apy. This second analysis provides 
additional information, such as the weak form factor, and 
clarifies the (modest) model assumptions necessary to ex- 
tract Rn- 

We start with distorted wave calculations of Any for an 
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electron moving in Coulomb and weak potentials f? . We 
use these to extract the weak form factor from the PREX 
measurement. In Born approximation, one can deter- 
mine the weak form factor directly from the measured 
Apy. However, Coulomb distortions may make Apy sen- 
sitive to the weak form factor for a range of momentum 
transfers q. In addition, the experimental acceptance for 
PREX includes a range of momentum transfers for labo- 
ratory scattering angles from about 3.5 to 8 degrees [5]. 
Therefore we will need to make very modest assumptions 
about the shape of the weak form factor (how it depends 
on momentum transfer q) in order to determine the value 
of the form factor at the average momentum transfer q 
0, 



{Q 



2x1/2 



0.475 ± 0.003 fm" 



(1) 



We initally assume the weak charge density of ^"^Pb, 
Pwif) has a Wood Saxon form. 



Pw{r) = 



1 



Po_ 

cxp[(r - 



R)/a] 



(2) 



TABLE I: Least squares fits of Wood Saxon {R, a, see Eq. [2| 
or Helm model {Ro, cr, see Eq. [sjl parameters to theoretical 
mean field model weak charge densities. 

Wood Saxon Helm 
Mean field force R (fm) a (fm) _Ro (fm) o" (fm) 

Skyrme I [S] 6.655 0.564 6.792 0.943 

Skyrme HI [9] 6.820 0.613 6.976 1.024 

Skyrme SLY4 [TOl 6.700 0.668 6.888 1.115 

FSUGold [TT] 6.800 0.618 6.961 1.028 

NL3 [H] 6.896 0.623 7.057 1.039 

NL3p06 [6] 6.730 0.606 6.886 1.010 

NL3m05 [6] 7.082 0.605 7.231 1.012 



Average 



0.61 ± 0.05 



1.02± 0.09 



neutron radii The average value of a for these mod- 
els is 0.61 ± 0.05 fm. Using a central value of a = 0.6 fm 
we obtain, 

Fwiq) = 0.204 ± 0.028(exp) ± O.OOl(mod). (6) 



with parameters po, R and a. Note, this form is only 
used to access the sensitivity to the shape of the form 
factor and our results will be independent of this assumed 
form. The weak density is normalized to the weak charge 
Qw = J (Prpwir), see below. 

We define the weak form factor F^^^q) as the Fourier 
transform of pw 



Fwiq) 



1 

Qw 



3 smgr 

d r pw[r)- 

qr 



(3) 



This is normalized Fwiq = 0) = 1. Our procedure is 
to calculate Ap^iO), including full Coulomb distortions 



[2], assuming pw from Eq. ^ 
laboratory scattering angle 0~i 
ceptance e(0) |5], 



We average ApyiO) over 
using the experimental ac- 



(A) 



J desmeeie)'i^A. 



J d9 sine ei9) 



(4) 



Here the unpolarized elastic cross section is ^ . We then 
adjust R until the calculated (A) agrees with the PREX 
result [3| 

A^^ = 0.656 ± 0.060(stat) ± 0.014(syst) ppm. (5) 

Here the first error is statistical and the second error 
includes systematic contributions. For a = 0.6 fm, we 
obtain a central value of i? = 6.982 fm, see below. Fi- 
nally from the pwii^) in Eq.j2l that reproduces Apy, we 
calculate F^iq) using Eq. ^ This procedure fully in- 
cludes Coulomb distortions and depends slightly on the 
assumed surface thickness a in Eq. [2] In Table |l] we show 
Wood Saxon fits to seven nonrelativistic and relativistic 
mean field model weak charge densities considered in ref. 
[B]. Note that these models span a very large range of 



Here the first experimental error is from adding the sta- 
tistical and systematic errors in Eq. [5]in quadrature. The 
second model error is from varying a by ±0.05 fm. This 
shows that the extracted form factor is all but indepen- 
dent of the assumed shape of the weak charge density. 
Equation [6] is a major result of this paper. This is the 
form factor of the weak charge density that is implied by 
the PREX measurement. 

We now explore some of the implications of Eq. |6] us- 
ing the Helm model [13) for the weak form factor. In the 
past, the Helm model has proven very useful for analyz- 
ing (unpolarized) electron scattering form factors [14. .15], 
see also ref. [16] for an application of the Helm model to 
neutron rich nuclei. The weak charge density is first as- 
sumed to be uniform out to a diffraction radius Rq. This 
uniform density is then folded with a gaussian of width a 
to get the final weak density. The width a includes con- 
tributions from both the surface thickness of the point 
nucleon densities and the single nucleon form factor. In 
the Helm model, the weak form factor has a very simple 
form, 



Fwiq) = —^jiiqRo)e 
qRo 



(7) 



with ji (x) — sin xjx^ — cos x/x & spherical Bessel func- 
tion. The diffraction radius i?o determines the location 
go of the zero in the weak form factor F-^^iq^) = 0. In 
coordinate space, the Helm model weak charge density 
can be written in terms of error functions (erf). 



pwir) = 



3Q 
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^(^0}. (8) 
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The root mean square radius of the weak charge density 
Rw (or weak radius) is i?^ = / d'^rr^ Pw{i')/Qw, 



Rw — 



[Ri + ha' 



(9) 



We see that Eq. [6] imphes via Eq. [7] a relationship be- 
tween allowed values of i?o and a. This relationship then 
implies via Eq. [9] a range of weak radii. Thus Eq. [6] does 
not, by itself, determine the weak radius. In principle 
the rms radius follows from the derivative of the form 
factor with respect to at g = 0. Because the PREX 
measurement is at finite g, one needs to assume some 
information about the surface thickness a in order to ex- 
tract Alternatively within the Helm model, if one 
determined the location of the zero of the form factor q^, 
in addition to Eq. |6j then this would uniquely fix both 
i?o and (7 and so determine Rw- 

In Table |T] we collect values of a determined by least 
squares fits of the Helm density, Eq. |8j to seven model 
mean field densities. The average of a for the seven mean 
field densities is 1.02 fm and individual results deviate by 
no more than 0.09 fm from this average. If one assumes 
cr = 1.02 ±0.09 fm, Eqs. [6) [7] and [9] imply 

Rw = 5.826 ± 0.181(exp) ± 0.027(mod) fm. (10) 

Again the larger experimental (exp) error is from adding 
the statistical and systematic errors in Eq. [5] in quadra- 
ture, while the model (mod) error comes from the coher- 
ent sum of the assumed ± 0.09 fm uncertainty in cr and 
the ±0.001 model error in Fw- The model error in Eq. 
10 provides an estimate of the uncertainty in Rw that 
arrises because of uncertainties in the surface thickness. 
Of course, it is not guaranteed that all theoretical models 
will have a surface thickness within the range 1.02 ± 0.09 
fm. Nevertheless, this result suggests that uncertainties 
in surface thickness are much less important for Rw than 
either the present PREX experimental error or even that 
of an improved measurement where the experimental er- 
ror is reduced by about a factor of three [17] . This is 
consistent with earlier results of Furnstahl suggest- 
ing a nearly unique relation between Fw (9) and the point 
neutron radius i?„ . We emphasize that if uncertainties in 
the surface thickness are a concern, one should compare 
theoretical predictions for the form factor Fwiq) to Eq. 
[6j instead of comparing theoretical predictions for Rw to 
Eq. [To] 

Comparing Eq. [TO] to the experimental charge radius 
Rch — 5.503 fm [21 [TU] implies a "weak charge skin" of 
thickness 



TABLE II: Helm model weak charge density parameters i?o 
and (J that reproduce the following values for the weak form 
factor Fwiq), see Eqs. [6]and[7] 



Density 


Ro (fm 


) cr (fm 


) Fw{q) 


Central value 


7.167 


1.02 


0.204 


Exp error bar 


7.417 


1.02 


0.176 


Exp error bar 


6.926 


1.02 


0.232 


Model error bar 


7.137 


1.11 


0.203 


Model error bar 


7.194 


0.93 


0.205 



that the weak charge density, of a heavy nucleus, is more 
extended than the electromagnetic charge density. 

In Fig. [l]we show a Helm model weak charge density 
that is consistent with the PREX measurement. This 
figure shows an uncertainty range from the experimental 
error and a model uncertainty from the assumed ±0.09 
fm uncertainty in a. Parameters for these densities are 
presented in Table [H] We also show in Fig. [ijthe (elec- 
tromagnetic) charge density [3 and a typical mean field 
weak charge density based on the FSUGold interaction. 



see Eq. 17 below. This theoretical density is within the 
error bars of the Helm model density. 
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FIG. 1: (Color on line) Helm model weak charge density 
~Pw{r) of ^"^Pb that is consistent with the PREX result 
(solid black line). The brown error band shows the incoher- 
ent sum of experimental and model errors. The red dashed 
curve is the experimental (electromagnetic) charge density p^h 
and the blue dotted curve shows a sample mean field result 
based on the FSUGold interaction [TT]. 



Rw - Rch = 0.323 ± 0.181(exp) ± 0.027(mod) fm. (11) 

Thus the surface region of ^°^Pb is relatively enhanced 
in weak charges compared to electromagnetic charges. 
This weak charge skin is closely related to the expected 
neutron skin, see below. Equation |11[ itself, represents 
an experimental milestone. We now have direct evidence 



Finally we wish to extract i?„ for ^"^Pb from Rw in Eq. 
[lOj We start by reviewing the relationship between the 
point proton radius Rp and the measured charge radius 
Rch- Ong et al. have [50] 

Rlh = Rl + (4) + ^{rl) + ^ + {r^)so- (12) 
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Here the charge radius of a single proton is {r"^) = 0.769 
fm^ and that of a neutron is (r^) = —0.116 fm'^. We 
calculate that the contribution of spin-orbit currents to 
Rch is small because of cancelations between protons 
and neutrons {r'^)so — —0.028 fm^. Finally the Darwin 
contribution 3/4M^ is also small with M the nucleon 
mass. For ^ospb we have, Rl,^ = i?p + 0.5956 fm^ or, for 
Rch = 5.503 fm [31112], 



Rp = 5.449 fm. 



(13) 



For the weak charge density of a spin zero nucleus, we 
neglect meson exchange and spin-orbit currents and write 
01 

pw{r)=^l dV[G,f(|r-r'|K(/) + G,^(|r-r'|)p,(r')] . 



(14) 

Here the density of weak charge in a single proton (r) 
or neutron (r) is the Fourier transform of the nucleon 
(Electric) Sachs form factors G^{Q^) and Gf (Q^). These 
describe the coupling of a boson to a proton or neu- 
tron m, 



4G^ = qp 



QnG^ — G 



(15) 



(16) 



At tree level, the weak nucleon charges are q^^ — —1 and 
Qp = 1 — 4sin^ Qw- We include radiative corrections by 
using the values (?„ = —0.9878 and qp = 0.0721 based 
on the up Gi„ and down Cid quark weak charges in refs. 
PSI [25] . The Fourier transform of the proton (neutron) 
electric form factor is G^(r) (G^(r)) and has total charge 
J (frG^Eir) = 1 (/ (frG%{r) = 0). Finally G% describes 
strange quark contributions to the nucleon's electric form 
factor pm24) . Note that there may be some small un- 
certainty regarding the dependence of the radiative 
corrections. This uncertainty could change R\, see be- 
low, by a very small amount of order (1 + g„)(rp). 

Equation [T4| can be rewritten by using a similar ex- 
pression for pch 



Pw{r) = qpPchir) + J d^r' [g„(G^p„ + C^pp) - G%pb] 

(17) 

with pi, — pn + Pp. The weak charge of -^^^Pb is 

Qw^ J d^rpwir) = Nq„ + Zqp = -118.55. (18) 
From Eq. [17] we relate the point neutron rms radius i?„, 

to i?M/, 



Rrn 



Qw n2 

it 

quN 



w 



quN 



Z 



Z + N 
qnN 



where {r1) = J d?r'r''^G%{r') is the square of the nucleon 
strangeness radius. This yields 

Rl = 0.9525i?^ - 1.671(r2) + 0.7450 fm^ (20) 
The strangeness radius of the nucleon (r^)^^^ is con- 
strained by experimental data [2TI - E1] and their global 
analysis [571 [55] • Using Table V of ref. [IB] for < q.II 
GeV2, gives (r^) = -6dG%/dQ^ = 0.02 ± 0.04 « ±0.04 
fm^. 



The neutron radius then follows from Eq. 10 



(19) 



R„ = 5.751 ± 0.175 (exp) ± 0.026(mod) ± 0.005(str) fm. 

(21) 

Here the very small third (str) error is from possible 
strange quark contributions. The neutron radius i?„ is 
slightly smaller than Rw because of the nucleon's size. 
Finally, the neutron skin thickness is 

R^-Rp = 0.302±0.175(exp)±0.026(mod)±0.005(str) fm. 

(22) 

This result agrees, within the model error, with the result 
of ref. ^, Rn-Rp = O-^^tolt fm. The smah difference 
between the present result and ref. [S] arrises because 
of small limitations of the Helm model in representing 
theoretical mean field densities. For example the Helm 
model does not have the correct expodential behavior at 
large distances. However, we have clarified how the ex- 
traction of the neutron radius depends upon assumptions 
on the weak skin thickness a and we provide an explicit 
model error for the uncertainty in _R„ — Rp because of 
uncertainties in a. 

We now summarize our results. In this paper we use 
distorted wave electron scattering calculations for 208pj^ 
to extract the weak charge form factor Fw{q), Eq. [6j the 
weak radius Rw, Eq- |10| and the point neutron radius 
Rn, Eq. [21] from the PREX parity violating asymme- 
try measurement. The weak form factor is the Fourier 
transform of the weak charge density at the average mo- 
mentum transfer of the experiment. This quantity is es- 
sentially model independent and is insensitive to assump- 
tions about the surface thickness. 

The extraction of Rw depends on modest assump- 
tions about the surface thickness. We use the Helm 
model to derive an estimate on the uncertainty in Rw 
because of the uncertainty in surface thickness. We find 
a "weak charge skin" where the surface region is rela- 
tively enriched in weak charges compared to electromag- 
netic charges. This is closely related to the neutron skin 
where i?„ is larger than the point proton radius Rp. Fi- 
nally, we extract given Rw, and find it to be slightly 
smaller than Rw because of the nucleon's size. 

We thank Witek Nazarewicz for very helpful discus- 
sions. We gratefully acknowledge the hospitality of 
the University of Tennessee and the Physics Division of 
ORNL where this work was started. This work was sup- 
ported in part by DOE grant DE-FG02-87ER40365. 



5 



[1] T. W. Donnelly, J. Dubach, Ingo Sick, Nucl. Phys. A503, 
589 (1989). 

[2] C. J. Horowitz, Phys. Rev. C 57 , 3430 (1998). 
[3] B. Frois et al., Phys. Rev. Lett. 38, 152 (1977). 
[4] C. J. Horowitz, S. J. Pollock, P. A. Souder, R. Michaels, 

Phys. Rev. C 63, 025501 (2001). 
[5] S. Abrahamyan et al. [PREX collaboration], Phys. Rev. 

Let. in press, 2012. 
[6] S. Ban, C. J. Horowitz, R. Michaels, J. Phys. G 39, 

015104 (2012). 

[7] X. Roca-Maza, M. Centelles, X. Vinas, M. Warda, Phys. 

Rev. Lett. 106, 252501 (2011). 
[8] D. Vautherin, D. M. Brink, Phys. Rev. C 5, 626 (1972). 
[9] M. Beiner, H. Flocard, N. Van Giai, P. Quentin, Nucl. 
Phys. A 238, 29 (1975). 
[10] E. Chabanat, P. Bonche, P. Haensel, J. Meyer, R. Scha- 

effer, Nucl. Phys. A 635, 231 (1998). 
[11] B.C. Todd-Rutel, J. Piekarewicz, Phys. Rev. Lett. 95, 
122501 (2005). 

[12] G. A. Lalazissis, J. Konig, and P. Ring, Phys. Rev. C55, 
540 (1997). 

[13] R. H. Helm, Phys. Rev. 104, 1466 (1956). 
[14] M. Rosen, R. Raphael, H. UberaU, Phys. Rev. 163, 927 
(1967). 

[15] R. Raphael, M. Rosen, Phys. Rev. C 1, 547 (1970). 

[16] S. Mizutori et al., Phys. Rev. C 61, 044326 (2000). 

[17] P. A. Souder et al., PREX H exper- 
imental proposal Jefferson Laboratory 
( http: / /hallaweb.jlab.org/parity/prex/prexILpdfj) . 

[18] R. J. Furnstahl, Nuc. Phys. A706, 85 (2002). 



[19] G. Fricke et al., 82-Pb Lead in H. Schoopper 
(ed.). Springer Materials -The Landolt- Born- 
stein Database ( http: / /www.springermaterials.com| ) 
DOLlO.1007/10856314-84. 

[20] A. Ong, J. C. Berengut, V. V. Flambaum, Phys. Rev. C 
82, 014320 (2010). 

[21] R. D. McKeown, Phys. Lett. B 219 (1989) 140.; D.T. 
Spayde, et.al. Phys. Lett. B583 (2004) 79; T. Ito, et.al. 
Phys. Rev. Lett. 92 (2004) 102003. 

[22] K.A. Aniol et.al, Phys. Lett. B 509 (2001) 211; K.A. 
Aniol et.al., Phys. Rev. C69 (2004) 065501; K.A. An- 
iol et.al, Phys. Rev. Lett. 96 (2006) 022003; K.A. An- 
iol et.al, Phys. Lett. B635 (2006) 275; A. Acha et.al, 
Phys. Rev. Lett. 98 (2007) 032301. Z. Ahmed et.al, 
'arXiv: 1107.0913 [nucl-ex]. 

[23] D.H. Beck, Phys. Rev. D 39 (1989) 3248; D.S. Armstrong 
et.al, Phys. Rev. Lett. 95 (2005) 092001; D. Androic 
et.al, Phys. Rev. Lett. 104 (2010) 012001. 

[24] F.E. Maas et.al, Phys. Rev. Lett. 93 (2004) 022002; F.E. 
Maas et.al, Phys. Rev. Lett. 94 (2005) 152001; S. Bau- 
nack et.al, Phys. Rev. Lett. 102 (2009) 151803. 

[25] J. Erler, A. Kurylov, M. J. Ramsey-Musolf, Phys. Rev. 
D 68, 016006 (2003). 

[26] K. Nakamura et al., JPG 37, 075021 (2010) 
(http://pdg.lbl.gov') Chapter 10 section 10.3. 

[27] R.D. Young, J. Roche, R.D. Carhni, A.W. Thomas, Phys. 
Rev. Lett 97 (2006) 102002. 

[28] J. Liu, R.D. McKeown, M.J. Ramsey-Musolf, Phys. Rev. 
C76 (2007) 025202. 



